Electron Delocalization in Gate- Tunable Gapless Silicene 
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Applying a perpendicular electric field can drive silicene into a gapless state, characterized by 
two nearly fully spin-polarized Dirac cones owing to both relatively large spin-orbital interactions 
and inversion symmetry breaking. Here we argue that since inter-valley scattering from non- 
magnetic impurities is highly suppressed by time reversal symmetry, the physics should be effectively 
single-Dirac-cone like. Through numerical calculations, we demonstrate that there is no significant 
backscattering from a single impurity that is non-magnetic and unit-cell uniform, indicating a sta- 
ble delocalized state. This conjecture is then further confirmed from a scaling of conductance for 
disordered systems using the same type of impurities. 

PACS numbers: 71.23.-k, 73.21.-b, 73.43.Nq, 
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It is well-known that single-cone Dirac fermion is im- 
mune to backscattering thus hard to be localized 043- 
However, graphene has two Dirac cones (valleys), as re- 
quired by the fermion doubling theorem^, H, 01 ■ Con- 
sequently, in the presence of impurities, the inter-valley 
scattering cannot be strictly prohibited and this leads 
to remarkable backscattering, resulting in localization 
m 2D [i-d^. This is essentially different from three- 
dimensional topological insulators (3DTIs), with just one 
Dirac cone for each surface [lll |. 

Recently, silicene, the silicon version of graphene on a 
honeycomb lattice, has been an exciting subject 12h15|. 
Due to its buckled structure, the spin-orbital coupling 
(SOC) is highly enhanced. With a perpendicular exter- 
nal electric field such structure also provides the tunabil- 
ity of the bulk gap Ag. As the applied field increases, 
the gap closing and reopening indicates a topological 
phase transition between a 2DTI and a trivial band in- 
" IM3. Exactly in the gapless critical state, 
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sulator 

where Aq = 0, the low-energy electronic structure can 
be described by massless Dirac Hamiltonian, forming two 
Dirac cones. The presence of various SOC interactions 
on the lattice results in rich spin textures around the 
Dirac points and eventually leads to profound behaviors 
in response to impurity scattering. 

The most intriguing property of the gapless gated sil- 
icene, also the focus in this work, is the opposite spin po- 
larization at different valleys, i.e., the valley-spin locking 
17 19]. Explicitly, the Dirac cone around K (K 1 ) point 
is polarized with spin up (down), mainly originating from 
the intrinsic SOC between next nearest-neighbor (NNN) 
sites as well as broken inversion symmetry (IS) due to 
external electric field. Thus, such pha se is dubbed spin- 
valley-polarization metal (SVPM) [18j. Ideally assuming 
no Rashba SOC, the spin around each cone is fully polar- 
ized, and, contrary to graphene, inter-valley (also spin- 
flip) scattering from non-magnetic impurities is strictly 



prohibited by time reversal symmetry (TRS). Therefore, 
two Dirac cones in this system are effectively decoupled 
and consequently the two-component, single- flavor Dirac 
physics emerges. Now it is quite essential to ask if there 
can be any delocalized states in strict sense under dis- 
order? In addition, Rashba SOC, which includes spin- 
flip processes, is nevertheless inevitable in realistic sil- 
icene. Can it induce inter-valley scattering and lead to 
the breakdown of the single Dirac cone physics as well? 

To answer these questions, in this paper, we systemat- 
ically study the non-magnetic impurity scattering prob- 
lem in a gapless system, designed to capture the physics 
of silicene and related materials, via numerical calcula- 
tions. By comparing with various typical arrangements 
of SOCs, we found that 1) from the quasi-particle in- 
terference (QPI) pattern associated with single impurity, 
within a certain region of parameter space (low energy, 
small Rashba SOC, and moderate impurity scattering 
strength) for SVPM, the "unit-cell impurity" will not 
give rise to significant inter- or intra- valley backscatter- 
ing; 2) the positive beta function (defined below) in the 
disordered system further confirms the conclusion in 1) 
and suggests the existence of a truly delocalized state. 

Hamiltonian and effective theory. Silicene or Ge, Sn, 
Pb counterparts can be minimally described by a four- 
band tight-binding model with energy scales t 3> Aso > 
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The first term describes the nearest-neighbor (NN) hop- 
ping with t = 1.6eV for silicene, where c\ a creates an 
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FIG. 1: (Color online) 2D dispersions E(k) (the upper row) and corresponding QPI curves 5p(Q) (the lower row) for IDF 
[(a) and (e)], graphene [(b) and (f)], and silicene with either Ar = [(c) and (g)] or Ar = Q.lt [(d) and (h)]. Red arrows on 
dispersions illustrate the orientations of pseudo-spins (a) or physical spins (the rest). Green thick arrows indicate significant 
scattering processes. Main QPI curves are plotted along the Q^-axis, while the inset of (g) in different directions. All QPI 
curves are plotted for impurity strength Vb = t, at Fermi energy Ep — 0.2t, with energy broadening 7 = 0.005t and 1000 x 1000 
grid for numerical integrations. 



electron at site i with spin polarization a. The second 
term represents the intrinsic SOC between NNN sites, 
where s = (s x , s y , s z ) are the Pauli matrices for physical 
spins, and v^j = (d% x dj) z /\di x dj\ = ±1 with and 
dj the two NN bonds connecting NNN sites i and j. The 
third term is the NNN Rashba SOC [13, where % = ±1 



for the A and B sites, respectively, and di 



dij j I d^ 



representing the unit vector of dij which connects NNN 
sites i and j. The fourth term represents the staggered 
potential, and the strength A v = l z E z can be tuned by 
a perpendicular electric field E z because of the buckling 
distance l z between two sublattices. If we take only the 
first term with t — 2.7eV, Eq. §T§ can also describe un- 
doped graphene. Hereafter, we adopt t as the energy unit 
and lattice constant a (NNN distance) as the length unit. 

Around two Dirac points at K(K') = (±47r/3, 0) in 
/c-space, the low-energy effective Hamiltonian for Eq. ([1]) 
with the basis (ipAf, ^Bf, i>Al, ~^b\) T reads 



(K(q) 9l (q) 

^t(-rjq x + iq y 

Ar^z + q y ) 






for pseudo-spin, with Fermi velocity ^t. Gating the 
system such that A v = Aso but with Ar, = 0, the full 
spin polarization of the valleys can be clearly seen: At 
valley K, the spin-up bands are gapless forming a Dirac 
cone, in contrast to spin-down bands now separated by a 
gap 2|Aso + A v | and thus out of the low-energy regime; 
at valley K' it is in opposite orientation due to TRS. 
The presence of considerable Ar destroys this full spin 
polarization but the majority around each valley would 
not change. Such states with two massless Dirac cones 
will be the main focus throughout this work. Restrict- 
ing A v = Aso while allowing to vary their strengths as 
well as the values of the Fermi level and Ar in the sys- 
tem give rise to rich physics, which reflects the interplay 
among spin, sublattice (pseudo-spin), and valley degrees 
of freedom under non-magnetic impurity scattering. 

QPI from single impurity. -We first investigate the 
scattering from a single impurity, by calculating QPI pat- 
tern [2l|. The Green's function for the clean system is 
G°(E,k) = G°{E,k,k) = [(E + - iJ(fc)] -1 , where 
I is the identity matrix and 7 <C 1 is the energy broad- 
ening. Here we only consider a single impurity with po- 
tential ~ 5(x) in a definite unit cell so that the impurity 
matrix V(k±, £2) = V 1S independent of k. The impurity 
induced Green's function is expressed as 

5G(E, k u k 2 ) = G°(E, kx)T{E, k 1 ,k 2 )G°(E, k 2 ). (5) 



where the wave vector q is measured from the Dirac 
point, rj = ±1 for K (K') point is the valley index, and 
h-iiq) is just the ideal Dirac fermion (IDF) Hamiltonian 



The standard perturbation method gives [2l| 

t(e) = [i -vr°(s)] _1 v; (6) 
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FIG. 2: (Color online) QPI pattern for silicene with different 
Fermi energies Ef, impurity strengths Vq, and Rashba SOC 
Ar. The upper row (a-c): Vo = t and Ar = 0; The middle 
row (d-f): Ef = 0.2t and Ar = 0; The lower row (g-i): Ef = 
0.2t and Vo = t. The scanning angle is chosen along the Q y 
axis from (a) to (f), where possible intra- valley backscattering 
reaches its maximum amplitude [3|], and along Q x axis from 
(g) to (i) for the detection of possible inter-valley scattering. 
All other parameters are the same with Fig. QJg). 
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G°(k,E). Now the Fourier trans- 



form of the induced local density of states is 

d 2 k 
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where Q = k' k and g(E,k,Q) = Tr(SG(E, fe, fc') - 
SG*{E,k',k)). The spectrum Sp(E,Q) in Eq. © is 
called the QPI pattern, which can also be obtained ex- 
perimentally from the Fourier transformation of STM 



measurements 
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This pattern provides an intu- 
itive picture of scattering processes: significant scattering 
processes will manifest themselves as peaks in the QPI 
pattern with associated scattering momenta. 

As a warm-up but essential example, we start with the 
single valley, single spin, 2x2 IDF Hamiltonian with just 
linear terms, /i^ -1 (fc). The impurity potential in fc-space 



can be V a = V t , V b = V t 3 , 
relative weight r, 



or their combinations with 



V = r ■ V a + (1 - r) • V b , 



< r < 1 



(8) 



with Ti the Pauli matrix acting on sublattice (pseudo- 
spin) space. The computed QPI, Sp(Q) of V a , is plotted 
as the blue curve in Fig. [He). The curve has no sig- 
nificant scattering peaks, consistent with the well-known 
fact that V a cannot induce backscattering for an IDF 0- 
,5J. Note that V a corresponds to a "unit-cell impurity" 



which is uniform within two sites of a unit cell. On the 
other hand, we also show the QPI for the "site impurity" , 
Vb, in Fig. 0Je) as the black curve. Two peaks associ- 
ated with intra-valley backscattering can be seen. This 
is not surprising because Vb is a mass term for IDF and 
destroys the pseudo- "TRS" , leading to a tendency to- 
wards localization^. From results for an impurity with 
different weights of V a and Vb also in Fig. [T^e) , it is in- 
teresting to notice that, a small weight (r > 20%) of 
V a is sufficient to annihilate the significant backscatter- 
ing peaks into the smooth background. In real space, 
an impurity with finite V a component corresponds to a 
long range one, with a smooth potential configuration 
within the unit cell. Such impurities can be dominant 
in graphenejH and therefore should also be easily real- 
ized experimentally for silicene. In the rest of the paper, 
we will restrict ourselves to the discussions of unit-cell 
impurity V a . 

Let us next consider the full tight-binding 4x4 Hamil- 
tonian H(k) in Eq. (J5J) with Aso = A v = Ar = 0, i.e., 
graphene without SOC. Compared with IDF Hamilto- 
nian, H(k) has two important features: Existence of 
two spins and two valleys, and higher order corrections 
(trigonal warping) within each valley, as illustrated in 
Fig. [IJb) ■ Given a unit-cell impurity potential, 



V a = V s <E> r = diag(V , V , V , Vo) , 



(9) 



which is "non-magnetic" both for physical spin s and 
pseudo-spin r, the corresponding QPI is shown in 
Fig. [Iff). It has very sharp peaks associated with inter- 
valley backscatterings between states with opposite k and 
velocity, as indicated by the green arrows. As in ordi- 
nary orthogonal disordered systems in 2D 25, |26|, this 
strong backscattering is responsible for the localization 
in graphene and weak 3DTI [lj]]. In short, the 

coupling between two Dirac cones (with opposite Berry 
curvatures [2?| ) makes the physics rather trivial. 

Armed with QPI patterns for above two examples, we 
come to our main focus, gapless silicene with Aso = 
A v = 0.5t. Such large SOC is taken simply for the pur- 
pose of giving enough space to extract out clear physics 
within our numerical precisions. No qualitative differ- 
ence is expected when inserting back the realistic value 
for silicene, as long as the Fermi energy lies within the 
window (— |Aso + A v |, |Aso + ^v|), where the SPVM pic- 
ture holds. We first consider Ar, = case. As mentioned 
before, this system is completely decoupled into two spin 
components and each valley is fully spin polarized [See 
Fig- IH c )]j an d therefore no inter- valley scatterings can 
happen. Indeed, given the same V a in Eq. (0), we notice 
that now the QPI in Fig. QJg) is qualitatively different 
from that of graphene, but similar to that of the IDF 
[Fig- [He)]. Moreover, there are also no significant intra- 
valley backscattering peaks. These qualitative features 
are isotropic in Q, as shown in the inset of Fig. [ljg), 
even though the full band structure is anisotropic around 
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each Dirac point. It has been argued that trigonal warp- 
ing would lead to nonzero backscattering amplitude [3|. 
However, our numerical results show that such backscat- 
tering is very weak and could be immersed in the contin- 
uum background of other scattering processes, reflected 
by the absence of a sharp and distinguishable peak in 
contrast to the case of graphene. Therefore, the gapless 
silicene effectively exhibits the massless Dirac fermion 
physics. Higher order corrections to IDF do not qual- 
itatively change the scattering behavior, so long as the 
Fermi energy is not far from the Dirac point and the 
impurity strength is not strong. This is one of the im- 
portant findings in this work. The absence of remarkable 
backscattering should signify a delocalized state to dis- 
orders, as will be numerically verified below. 

Before entering into the discussion on disordered sys- 
tems, two remarks are in order. First, in Figs. [TJc) 
and (g), with vanishing Ar, inter- valley scattering is in 
fact suppressed from a priori. Nonzero Ar, as to be the 
case in silicene, should couple states with opposite spin 
and makes the spin-valley polarization imperfect [See 
Fig. (Ud)]. However, as shown in Fig. [ljh), it is remark- 
able to see that such Rashba term does not give rise to a 
significant inter-valley scattering, and thus the effective 
"single-valley Dirac physics" remains intact. Second, in- 
creasing parameter such as Ef , impurity strength Vq, or 
Ar in the system is expected to enhance intra- and inter- 
valley scattering processes due to unavoidable contribu- 
tions from higher order corrections and spin/ valley mix- 
ing. Indeed, as clearly shown in Fig. [2j the QPI pattern 
changes at some point, indicating a transition from a de- 
localized to localized state beyond effective single- valley 
Dirac physics. For instance, in the case of very strong 
Ar, in Fig. [2] (i), although two states |fc) and | — fc) in 
different valleys (with exactly opposite spin orientations) 
cannot be coupled by a non-magnetic impurity, but spin 
orientations in their neighborhoods will not be exactly 
opposite. Thus an inter-valley backscattering can be al- 
lowed due to the energy broadening 7. 

Scaling of conductance: multiple impurities. So far, 
the scattering from a single impurity has been investi- 
gated. If the backscattering is effectively ignorable, does 
this really lead to delocalized ground state in disordered 
gapless silicene with unit-cell impurities? To confirm 
this, we perform a standard numerical scaling for dis- 
ordered silicene. Disorder is added to the Hamiltonian 
Eq. ([T]) as J2i a e i c ia c icri where ti is a random number 
uniformly distributed in (— W/2, W/2). Here q is inde- 
pendent of spin due to TRS. If for two sites in each unit 
cell are identical but randomly assigned values among 
different unit cells, then it corresponds to unit cell im- 
purities V a . The intrinsic conductance g is defined as 
l/g = 1/gL — l/-^c) where g^ is the two-terminal quan- 
tum conductance, N c is the number of propagating chan- 
nels and 1/N C is the contact resistance 28]. This g is 



suitable for a numerical scaling 



rfflng) 
dlnL ' 



(10) 



where (• • • ) is the average over random ensemble, and 
L is the spacial size of the sample with a fixed ratio of 
length and width. This scaling function f3 is used as a 
criteria: (3 < and /? > correspond to localized and 
delocalized states, respectively. 

In Fig. [3l we plot (In g) as a function of size L (in log- 
arithmic scale), where the slope represents f3. It can be 
seen from Figs. [3Ja) and (b) that, for unit cell impuri- 
ties, apart from some fluctuations due to the smallness 
of conducting channels, (ln<?) is clearly increasing with 
increasing L, suggesting delocalized state with (3 > 0. 
These are consistent with our results of the absence of 
significant backscattering from the single impurity study, 
further confirming the robustness of the effective single- 
valley Dirac physics. Note that this is totally different 
from the case of graphene with the W = 2t [black curve 
in Fig. HJc)], where the slope is negative. It has been 
found that for graphene, even long-range impurities can- 
not maintain a fully delocalized state with /3 > be- 
cause of inevitable inter- valley scattering • Of course, 
as in any lattice models, sufficiently strong disorder will 
eventually localize all the electrons, as the red curve in 
Fig. 12c) shows. Therefore, it is natural to expect rich 
localization-delocalization transition behavior in the pa- 
rameter space spanned by Ep, W ', Ar, and Aso (= A v ). 
More details of such localization-delocalization transi- 
tion, e.g., the universality, critical exponents, and global 
phase diagram will be discussed elsewhere. 

Delocalized bulk states in gapped Kane-Mele model 
with nontrivial Z2 topological nature were found in 
Ref. [24j ■ In that case, delocalized states can only appear 
when NN Rashba SOC is nonzero, i.e., when the system 
becomes symplectic; otherwise, the system is decoupled 
into two gapped unitary subsystems, where no states with 
/3 > can be observed 24j. This is indeed reasonable as 
a gapped Dirac cone has serious backscattering Q. In 
our case, however, the physics behind derealization lies 
on either independent unitary subsystems, each of which 
owns a gapless Dirac cone (zero Rashba SOC), or sym- 
plectic subsystems with two nearly independent gapless 
Dirac cones (nonzero Rashba SOC). 

In summary, we reveal the essential transport proper- 
ties via numerical simulations on a critically gated buck- 
led honeycomb structure of silicene (and also suitable for 
Ge, Sn, Pb counterparts) under non-magnetic impurity 
scattering. QPI by a single unit-cell impurity shows no 
significant backscattering, suggesting an effective single- 
valley Dirac physics, in spite of weak trigonal warping. 
The robustness of such delocalized state is further con- 
firmed by the positiveness of the /3 function for a disor- 
dered system, even in the presence of Rashba SOC. This 
feature sheds a new light on constructing high mobility 
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FIG. 3: (Color online) Typical conductance as a function of the system size (in logarithmic scale), (a) Aso = A v = 0.5t, Ar = 0; 
(b) Ar = O.lt; (c) Two examples of localization: graphene with W = 2t (black), and silicene with Aso = A v = 0.5i, Ar, = 0, 
under strong disorder W — 4t (red). Each dot is the average over 2000 disorder samples with unit cell impurities. 



silicene-based electronic devices. We believe our result 
is also insightful to relevant systems such as a 2D M0S2 
[30| and a cold-atom system with arranged SOC[3lj|. 
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